Spectral Correlation in Incommensurate Multi- Walled Carbon Nanotubes 



m 
O 
O 

(N 
C 

a 

Hi 

in 



K.-H. Aim, 1 Yong-Hyun Kim, 2 J. Wiersig, 3 and K. J. Chang 2 ' 4 

1 School of Physics, Seoul National University, Seoul 151-747, Korea 
2 Department of Physics, Korea Advanced Institute of Science and Technology, Daejeon 305-701, Korea 
3 Max Planck Institute for Physics of Complex Systems, 
Nothnitzer Strasse 38, Dresden 01187, Germany 
* School of Physics, Korea Institute for Advanced Study, Seoul 130-012, Korea 

(Dated: February 1, 2008) 

We investigate the energy spectra of clean incommensurate double-walled carbon nanotubes, and 
find that the overall spectral properties are described by the critical statistics similar to that known 
in the Anderson metal-insulator transition. In the energy spectra, there exist three different regimes 
characterized by Wigner-Dyson, Poisson, and semi-Poisson distributions. This feature implies that 
the electron transport in incommensurate multi-walled nanotubes can be either diffusive, ballistic, 
or intermediate between them, depending on the position of the Fermi energy. 



-a 
c 
o 

(N 
> 

oo 

(N 
OO 

o 

<N 
O 
-i— > 



i 

C 

o 
o 



X 



PACS numbers: 03.65. 



72.15.-v, 73.23.-b 



Carbon nanotubes have attracted great attention due 
to their remarkable electrical and mechanical proper- 
ties Q. Many theoretical and experimental studies have 
demonstrated that single- walled nanotubes (SWNTs) ex- 
hibit ballistic electron conduction 0, and single- 
molecule devices utilizing semiconducting SWNTs have 
been realized 0. On the other hand, despite useful ap- 
plications of multi- walled nanotubes (MWNTs) []|, the 
transport properties of MWNTs are not well understood 
and even controversial; conductance measurements usin 
scanning probe microscope showed ballistic behavior 
while magnetoresistances measured for MWNTs on top 
of metallic gate indicated diffusive conduction 0, ■ 

In general, MWNTs have very complex electronic 
structure, so that a direct use of their transport prop- 
erties is severely hindered. If concentric carbon shells 
are especially incommensurate, the Bloch theorem is no 
longer valid, and the Landauer formula for conductance 
is not applicable because it is difficult to count the num- 
ber of conducting channels 0. In fact, several experi- 
ments on MWNTs indicated that carbon shells have of- 
ten different periodicities 0, However, electronic 
structure calculations have been mostly focused for com- 
mensurate multi-walled nanotubes so far. Here, instead 
of directly calculating conductances, we investigate the 
spectral properties of energy levels in incommensurate 
MWNTs. 

The spectral analysis of energy levels has proven to be 
a useful tool to probe the nature of eigenstates in disor- 
dered conductors [l2T |. In the diffusive metallic regime, 
the statistics of spectral fluctuations is described by the 
Gaussian orthogonal ensemble (GOE) of random-matrix 
theory In this regime, the distribution P(s) of en- 
ergy spacings between nearest levels is well fitted by the 
Wigner-Dyson surmise -Pgoe(s) = §s exp(— f s 2 ), where 
s is in units of mean level spacing A. In the insulating 
regime, where the energy levels are uncorrelated, P(s) 
is given by the Poisson distribution, Pp{s) — exp(— s). 



At the Anderson transition point, a critical statistics 
was found The short-range level correlation has 

been later described by semi-Poisson distribution [l2| . 
which has now been an important issue in quantum chaos 
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In this Letter we make a spectral analysis of the en- 
ergy levels of clean incommensurate MWNTs, and find 
that the level statistics is similar to that of the Ander- 
son metal-insulator transition. The energy spectra ex- 
hibit both the Wigner-Dyson and Poisson distributions 
depending on the position of the Fermi level. In ad- 
dition, we demonstrate that incommensurate MWNTs 
have an intermediate regime, which satisfies the semi- 
Poisson statistics. 

To investigate the spectral correlation, we consider 
double- walled carbon nanotubes (DWNTs), because the 
electrical conduction in MWNTs is believed to be gov- 
erned by the outermost shells 0, 0. We use a tight- 
binding model with one 7r-orbital per carbon atom, which 
successfully describes the electronic structure of DWNTs 
mm The tight-binding Hamiltonian is given by 
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where 70 (= -2.75 eV) is the hopping parameter between 
intra-layer nearest neighbor sites, i and j, and W (= 
70/8) is the strength of inter- wall interactions between 
inter-layer sites, i' and j', with the distance of d^f and 
the cut-off for > 3.9 A. Here 9ij is the angle between 
two 7T orbitals, <v is the annihilation operator of an elec- 
tron on site i', a (= 3.34 A) is the distance between two 
carbon walls, and S = 0.45 A. 

To extract the fluctuations from the level sequence, it 
is customary to map the real energy spectra {q} onto 
the unfolded spectra {Ei} through Ei = N(ci), where 
N(ei) is the number of levels u p to e, and the over- 
line denotes its broadened value |2lJ. Here we use the 
Gaussian broadening scheme for averaging the energies 
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FIG. 1: The nearest energy spacing distribution P(s) for 
energy levels between -7 and 7eV and different helicities 
(16,5)/(23,8) (triangle), (17,2)/(16,15) (circle), (16,5)/(17,15) 
(square), and (17,2)/(21,9) (diamond). The semi-Poisson 
(dashed), Poisson (long-dashed), and GOE (dotted) distribu- 
tions were also plotted for comparison. The inset shows the 
number variance S(L) for (16,5)/(23,8) case. The number of 
considered levels for (16,5)/(23,8) case is 20220. 



|22|. After unfolding, we obtain the distribution P(s) 
as a function of nearest energy spacing, Si = Ei — E i _ 1 
(the mean level spacing of the unfolded spectra equals 
1), and the number variance X(i-) = ((n(L,_E) — £) 2 ), 
which represents the variance of the number of levels 
in the interval (E - L/2,E + L/2), i.e., n(L,E) = 
N(E+L/2)-N(E-L/2). The bracket ( • ■ • ) denotes the 
average around the energy E on an energy window much 
larger than the mean level spacing but much smaller than 
E. 

The two shells of DWNTs are incommensurate if the 
ratio of the unit cell lengths along the tube axis is irra- 
tional. We test many incommensurate DWNTs, where 
both the shells arc metallic or semiconducting, a semi- 
conducting shell is inside a metallic tube, and vice versa. 
The diameters of the inner and outer tubes are 1.5 and 
2.2 nm, respectively, and the nanotube length is set to 
be about 53 nm. We find the spectral properties are sim- 
ilar for all the incommensurate DWNTs with different 
helicities. (See the P(s) in Figure P). Since the spec- 
tral properties are similar for all the tubes considered 
here, from here on, we only present the spectral proper- 
ties of the (16,5)/(23,8) DWNT where the semiconduct- 
ing (16,5) single-walled tube is aligned inside the metallic 
(23,8) tube. In the inset of Figure £(£) is also shown 
for the (16,5)/(23,8) nanotube. We find that P(s) and 
£(£*) cannot be described by the Poisson distribution or 
GOE, but, they are well described by the semi-Poisson 



(SP) distribution [3 HIS; 

P SP {s) = 4sexp(-2 S ) (2) 
Zsp(L) = | + i(l- e - 4L ). (3) 

The SP distribution is defined by removing every other 
level from an ordered Poisson sequence and turned out 
to be a reference point for the critical statistics of several 
disordered systems like the Anderson model at mobility 
edge 0, H3| , and also of other systems without disorder 
such as pseudo-integrable quantum billiards 0, 0, 0] . 

While P(s) carries information on short-range corre- 
lation in the energy spectra, contains rather long- 
range correlation. Usually X(£) indicates a deviation 
from the universal value for large L, thus, the particle 
dynamics becomes nonuniversal at short time scale. For 
large L, the spectral correlation is not universally semi- 
Poisson but depends on the helicity of nanotubes. The 
deviation from the SP distribution gives useful informa- 
tion; the linear behavior for large L, i.e., T,(L)/L — » x, 
represents the level compressibility, and for disordered 
metals \ was shown to be [25| 




where Di is the multifractal exponent of the inverse par- 
ticipation ratio and d is the spatial dimension. Here D 2 /d 
is related to the 'probability of return' [2(|, 

p(t) oc t-°^ d . (5) 

Our incommensurate DWNT behave as a two- 
dimensional systems (d = 2) with off-diagonal disorder, 
which has a mixed boundary between periodic and hard- 
wall conditions. From the inset in Fig. ^ we estimate 
D 2 /d to be about 0.32 for the DWNT considered here. 
While p(t) oc t~ d / 2 in normal diffusion, our system ex- 
hibits anomalous diffusion, which was also found by pre- 
vious wavepacket spreading calculations |2o| . 

The SP distribution in Fig. is mostly contributed 
from the region of large \E\, where the density of states 
(DOS) is high, while the energy statistics is usually not 
independent of the energy regime. Since real electron 
conductions occur near the Fermi level, it is instructive 
to examine the energy-level statistics on various energy 
windows. Although we discuss the level statistics when 
the Fermi level increases, we find similar statistics for the 
downward shift of the Fermi level, which usually occurs 
in hole-doped tubes. 

In the level statistics, we remove the states lying be- 
tween and 0.08 eV, which are associated with the lo- 
calized states near the sample boundary due to the finite 
size of the system. In Fig. EJa), P(s) is drawn for low en- 
ergy states from 0.08 to 1.3 eV, and exhibits the Poisson 
distribution, indicating that the energy levels are uncor- 
rected. Single-wall nanotubes, which are constituents 
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of the DWNT, are approximately described in terms of 
Pp(s) and S(s). In this case, the delta-function peak re- 
sults from the doubly degenerate states due to inversion 
symmetry. If inter-wall interactions are absent in the 
DWNT, a superposition of two independent spectra of 
the nanotube constituents is also described by the Pois- 
son and delta functions. 

When inter-wall interactions break the degeneracy, 
smearing out the delta-function peak, P(s) is still 
Poisson-like, as shown in Fig. 2(a), implying that the 
energy mixing between the two carbon shells is insignif- 
icant for low energies. In the low energy regime, it is 
useful to use the Landauer formula, 

G =X^ Ti( ^ } ' (6) 

i 

where Ti(Ep) is the transmission probability of the chan- 
nel i obtained from the band structure calculations at the 
Fermi energy Ep Our results infer that quantized 
conductances observed in MWNTs [|| may be due to the 
fact that the Fermi level is close to the charge neutrality 
point, Ep ks 0. However, it is usually very difficult to 
find experimentally the location of the Fermi level. 

As the Fermi level is shifted further, since the mean 
level spacing becomes smaller and the inter-wall interac- 
tion matrix elements become more significant, the inter- 
wall interactions induce a more significant mixing of the 
energy levels between the carbon shells. For higher en- 
ergy states, we find that P(s) follows the SP distribution, 
as shown in Figs. Hfb) and (d). The quantum conduc- 
tance in this regime still remains a challenging problem. 
In particular, it is a nontrivial but important task to find 
the tube length dependence of quantum conductance in 
this regime. Recent calculations [27} showed negligible 
inter- wall tunneling currents for the tube length of 1 //m, 
but it is questionable whether one can use the Bardeen's 
formalism for DWNTs where the coherent back-tunneling 
between carbon shells is non-negligible. According to the 
weak localization formula, 5a oc — J dtp(t), with Eq. J^J, 
conductances may depend on the tube length. However, 
the weak localization formula considers only the inter- 
ference of the pair of time-reversal paths, which is not 
guaranteed in the critical regime. 

In Fig. |2Jc), one can see that the critical statistics 
evolves to the GOE as going to the higher energy re- 
gion, where P(s) is quite close to the Wigner-Dyson sur- 
mise. We also find the GOE-like statistics for both P(s) 
and S(-L) in all the incommensurate DWNTs tested here. 
The appearance of the GOE-like behavior indicates that 
all the symmetries except for time-reversal symmetry are 
effectively broken, and the corresponding electron motion 
is ergodic. In this regime, the conductivity of MWNTs 
may be described by the weak localization theory. One 
should also note that the weak localization theory is ap- 
plicable for energies higher than those in the first critical 



Fig. 2 Ahn, et al. 
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FIG. 2: The nearest energy spacing distribution P(s) on 
the energy windows of (a) (0.08,1.3eV) with 511 levels, (b) 
(1.0,2.5eV) with 1872 levels, (c) (2.5,3.5eV) with 2696 levels, 
and (d) (3.5,4.0eV) with 888 levels. 



regime, because the sequence of statistics with increasing 
energy is Poisson — * SP — > GOE — > SP . 

In MWNTs on top of metallic gates, the details of 
contacts, local excessive charges, and depletion of charge 
carriers may change the position of the Fermi level |28| . 
Based on our results, we guess that in conductance mea- 
surements 0, S Hi , the Fermi level is shifted to the GOE 
regime or to the critical regime (close to the GOE) where 
conductivity is properly described by the weak localiza- 
tion formula. 

The calculated variances over every 1000 consecutive 
levels are compared with the densities of states for two 
different DWNTs in Fig. |3| The size of fluctuations of 
the level spacing is dictated by the spectral statistics; 

poo 

(( s -( s )f)= / ds(s-l) 2 P(s) (7) 
Jo 

= 1, i - - 1 (Poisson, SP, GOE). 

2 7T 

It is expected that the effect of inter-wall interactions on 
the level spacing is most significant when the mean level 
spacing A is minimum near E = ±70 = ±2.75eV, where 
the maximum DOS occurs on a graphene sheet. One can 
see ((s-(s)) 2 ) = |-1 « 0.27 near E= ± 7o = ±2.7beV, 
which is a signature of GOE distribution. The variances 
larger than 1 near E = are due to the fact that P(s) is 
not completely Poisson-like because the degenerate levels 
of each nanotube layer are not completely lifted up. 

We find that a transition from the Poisson-like to SP- 
like statistics occurs at lower energies as the diameter of 
nanotubes increases. One of the reasons for this trend 
seems that the DOS is enhanced for nanotubes with 
larger diameters. The transition energy is close to 1 eV 



4 



Fig. 3 Ahn, et al. 
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FIG. 3: The variance of the level spacing (solid) in Eq. 
and the density of states (dotted) in arbitrary units for the 
(16,5)/(23,8) (upper panel) and (43,8)/(51,9) (lower panel) 
double-walled tubes. 



in the (43,8)/(51,9) DWNT with the outer and inner di- 
ameters of 4.4 and 3.7 nm, respectively, while it is higher 
than 1 eV for the (16,5)/(23,8) tube with the diameters of 
2.2 and 1.5 nm, as shown in Fig. Since experimentally 
measured diameters of MWNTs are often in the range of 
10 nm, we expect that the crossover of spectral statistics 
occurs at energies lower than 1 eV. 

Very recently, Kodak and his co-workers |2^ measured 
the conductance of a DWNT whose two tubes have a gap 
and showed linear conductance near the Fermi level. The 
linear conductance or a finite DOS near the Fermi level 
indicates that the Fermi level indeed can be shifted far 
from the charge neutral point. 

Finally, we point out that the spectral statistics might 
be probed through conductance measurements in the 
Coulomb blockade regime[30l l3l| . The temperature 
should be low enough to ensure that ksT < U, A, where 
U is the charging energy of nanotube dot. In the resonant 
tunneling regime, where electron tunnels through one 
quantum level, if electron-electron interactions are not 
too strong 32], the conductance peak spacing with vary- 
ing the gate voltage provides information on the single- 
particle energy spacing. 

In conclusion, we have shown that the spectral statis- 
tics of incommensurate double-walled nanotubes follows 
the Poisson, GOE, or SP distribution, depending on the 



energy window, while the overall states are well described 
by the SP distribution. This results indicate that the na- 
ture of electron transport in multi-walled nanotubes can 
be either ballistic, diffusive, or critical, depending on the 
position of the Fermi level. It is questioned whether the 
usual weak localization correction is relevant to existing 
experiments. The Coulomb blockade oscillation in nan- 
otube dots is suggested to investigate the spectral statis- 
tics in this work. 
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